Quantum Particles Constrained on Cylindrical Surfaces with Non-constant Diameter 
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We present a theoretical formulation of the one-electron problem constrained on the surface of 
a cylindrical tubule with varying diameter. Because of the cylindrical symmetry, we may reduce 
the problem to a one- dimensional equation for each angular momentum quantum number m along 
the cylindrical axis. The geometrical properties of the surface determine the electronic structures 
through the geometry dependent term in the equation. Magnetic fields parallel to the axis can 
readily be incorporated. Our formulation is applied to simple examples such as the catenoid and 
the sinusoidal tubules. The existence of bound states as well as the band structures, which are 
induced geometrically, for these surfaces are shown. To show that the electronic structures can 
be altered significantly by applying a magnetic field, Aharonov-Bohm effects in these examples are 
demonstrated. 
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I. INTRODUCTION 

Developments in the synthesis of nano-structured ma- 
terials is ever expanding the list of structures which may 
be used for applications. Among other things, structures 
obtained as self-assembled interfaces in lyotropic liquid 
crystals, which are used as the templates of meso-porous 
materials, as well as structures made of bent graphites, 
including carbon nano-fuUerenes, are intriguing exam- 
ples. These structures are characterized by curved sur- 
faces in nanometer scales and often possess a variety of 
geometrical properties. It is interesting to study the elec- 
tronic properties related to the geometrical properties if 
the electron is constrained on the curved surfaces, be- 
cause such unique settings may provide a new phenomena 
in quantum domain which can be pursued for electronic 
applications. We here consider the problem of quantum 
particle motion constrained on curved surfaces with the 
tubular topology and cylindrical symmetry. 

We start with a model approach developed in Refs.jlj, 
0, Q) in which the electron is assumed to be confined 
strongly to within a thin layer built over the surface. A 
limiting procedure that brings the thickness of the layer 
to be zero leads to the Hamiltonian in the form of the 
sum of a Laplace-Beltrami operator and an effective at- 
tractive potential energy associated with the local curva- 
tures. These terms are written using double curvilinear 
coordinates on the surface. The stationary solutions of 
the Schrodinger equation are inevitably affected by the 
metric properties, curvatures, and topology of the sur- 
face. Previous studies reported that the effective poten- 
tial energy can induce bound states whose localization 
centers are strongly curved parts of the surface. 0, Q 
The discrete energy levels for bound states usually lie 
below the lowest edge of the essential (continuous) spec- 
trum. The electronic structures can also be affected sub- 
stantially by the global topology of the surface, because 
it restricts the ways in which the electronic wave propa- 
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gates and interferes. This point, however, has not been 
studied comprehensively. 

The purpose of this paper is to develop a basic formu- 
lation for curved surfaces with the tubular topology and 
cylindrical symmetry, and also to present some of its out- 
comes for simple examples. Our attention is focused on 
the case of tubular surfaces with non-constant diameter, 
namely, ones with non-vanishing Gaussian curvatures. 
The geometrical character of the surfaces dealt with in 
this paper is among the simplest and can be treated with 
relatively simple mathematics, though the results are not 
trivial. Within our formulation, it can be readily shown 
that there can be bound states lying inside the essential 
spectrum with a different angular momentum. The anal- 
yses of the electronic structures of curved surfaces with 
more complicated geometries, especially ones so called 
triply periodic minimal surfaces, should resort to more 
elaborated group theoretical techniques. Our formula- 
tion is expected to provide some insight into further stud- 
ies. 

This paper is organized as follows. In Section 2, the 
basic Schrodinger equation for an electron strongly con- 
strained on a curved surface is summarized briefly, l,*?,*^ 
The equation is applied to a class of tubular surfaces 
with cylindrical symmetry in Section 3. By symmetry, 
the basic two-dimensional equation can be reduced to a 
one-dimensional form, where an axial magnetic field can 
also be taken into account. In Section 4 we study several 
specific examples and provide the functional form of the 
potential energy for each case. These specific problems 
are solved numerically, where we observe the existence 
of bound states at the strongly curved region on a sur- 
face. For periodically modulated structures, the relevant 
one-dimensional band structures as well as charge den- 
sity distributions along the tubule are calculated. The 
electronic structures can be manipulated by tuning the 
applied magnetic field. The Aharonov-Bohm oscillations 
of the energy levels of the bound states arc demonstrated. 
Further discussions are given in Section 5, and the con- 
clusions in Section 6. 
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II. ONE-PARTICLE MOTION CONSTRAINED 
ON A CURVED SURFACE 

Consider a two-dimensional manifold Ai embedded in 
the three-dimensional Euclidean space and assume 
a curvilinear surface coordinate system (i — 1,2) on 
M. Then the metric tensor gij = 1,2) is defined by 
ds^ = gijdq^dq^ where ds is an infinitesimal distance. We 
write the principal curvatures as (i — 1,2). The quan- 
tum mechanics of a single electron which is constrained 
physically on M through a strong confining potential en- 
ergy is described by the following Schrodinger equation, 

1 



) = Ei;iq\q% (1) 



where g — det (gij ) and (g* 



. The first term 
involves the Laplace-Beltrami operator which describes 
the propagation of waves through a curved surface with 
non-Euclidean metric. The second term is the attractive 
potential energy associated with the local curvature, |5| 
which essentially arises from the 'spread' of the wave 
function perpendicular to the surface and thus is related 
to the uncertainty principle. The normalization condi- 
tion for ipil^ J Q^) is given by the integral 



mq\q^)\^gidq'dq^ 



(2) 



which should be unity for square integrable cases. Static 
magnetic fields can also be taken into account, in which 
case the equation becomes 



1 1 / a e 

~7.h 1- - , 



, 'ih-^ + -A^ ] .fqq''^ ( ~ih-^ + -A, 
2M^\ dq' c 7 \ dqi c ^ 



Xi^{q\q') - - K^fi'iq^q^) = E^{q\q^), (3) 

where Ai are the components of the magnetic vector po- 
tential associated with the coordinates g*. 



III. SURFACES WITH CYLINDRICAL 
SYMMETRY 

Now let us consider a class of smoothly curved surfaces 
with cylindrical symmetry, so that they are invariant un- 
der an arbitrary rotation around the z axis. It is con- 
venient to use the cylindrical coordinates (r, 0, z), with 
which the generator for a rotation around the z-axis is 
given by Iz = —id/d(j), the z-component of the angular 
momentum. Any of such surfaces can be specified by the 
set of coordinates {z{x),r{x)) as functions of the longi- 
tudinal surface coordinate x taken perpendicular to the 
azimuthal coordinate, (/>. A given value of x thus speci- 
fies a circle of radius r{x) lying on z = z{x) plane. If the 



fimction z{x) has a monotonic dependence so that it is 
invertible, then the radius can be expressed as a function 
of z, r = r(x{z)) — r{z). This means that, in describing 
such a surface, we may use the z coordinate as the longi- 
tudinal surface coordinate. We shall treat this particular 
case in more detail in the following. 

For the cylindrical surface specified as r = r(z), the 



surface metric ds^ = gzzdz'^ + 

9zz = 1 + r'izf 
g<k^ = r{zf, 



is given by 



(4) 
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r{zf{l + r'{zf). 



The principal curvatures of the surface are given by 
1 r"{z) 



K2 = 



r{z)Jl + r'{zy 



(l + r'(z)')i' 



(5) 



where K2 corresponds to the curvature along the longitu- 
dinal direction. 

We shall consider a magnetic field which is applied to 
the system parallel to the z-axis, and suppose its magni- 
tude is given by B(r), a function of the radial coordinate 
r. In the symmetric gauge, the vector potential is given 

by 



{Ax, Ay, Az 



$(r) 1 
2tt 



(6) 



where $(r) — 27rr' B{r')dr' is the total flux passing 
through a disc of radius r. In the cylindrical coordinates 
(r, (j), z), the relevant components of the vector potential 
are given by 



{Ar , Acf, , Az 



(0,^,0). 

ZTT 



(7) 



This can be applied for an Aharonov-Bohm type mag- 
netic field, A^ = $/27r, as well as a uniform magnetic 
field, A^ = Br'^/2. The Schrodinger equation can now 
be written as 



2M 



1 



d 



d 



12 dz 



XlP{z,(l3) 



K2fi^{z,4,)^Ei;{z,ct>),{S) 



where Iz = —i-§^- Since Iz commutes with the Hamil- 
tonian (the rotational symmetry), we can assume as the 
form of the solution 'ip{z, 0) = e""''^4'm{4')i in which m is 
the eigenvalue of Iz and can take only integer values in 
order for the wave function to be single valued. Then the 
Schrodinger equation is written as 
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Vl + 7-'^ dz y^ip: 
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(k1 - K2flpra{z) = Elprniz) , (9) 
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where 



IV. SPECIFIC EXAMPLES 



A™(z) = TO + —A^{r{z)). 
This has a Sturm-Liouvillc form, 



(10) 



d 



p{z)- 



dz 



+ q{z)^miz) = ew{z)'ipmiz), (11) 



where 



and 
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(12) 



piz) 
Qiz) 



w{z), 



(13) 
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w(2:) = rVl + r' 
Hence by introducing the new variable x as 



(14) 



Eq.0 can be reduced to a simpler form,[ 

d2 



da; 2 



y„(a;) + W™(x)j/™(x) = ej/m(x) 



where 

ym(a;) = {p{z)w{z))iijmiz) = ^r(z)'(/'m(^), 

and 



(15) 



(16) 



1 



g(^) , 

w{z) {p{z)w{z))i dx^ 
1\ 1 1 



X^mizY 



A J r(z)2 



{p{z)w{z))-^ (17) 
(18) 



Ea. (|15|l has the form of a one-dimensional Schrodinger 
equation with a potential energy Wm(cc), where x cor- 
responds to the Euclidean line length along a curve on 
the surface with fixed cj). Note that the second term in 
Eg. lfTTIl turns out to be (k^ — 2kiK2 — l/r^)/4, so that 
all terms including ki are formally cancelled out in the 
last expression (|18|) . The normalization condition for the 
wave function ym{x) is 



27r / \y^{x)\^dx^l, 



(19) 



for square integrable cases. 



We shall study several particular examples of curved 
surfaces with cylindrical symmetry by using the formula- 
tion presented above. The first and the simplest example 
is the catenoid, which gives one of the simplest examples 
of curved surfaces with zero mean curvature or minimal 
surfaces. The second example is the tubular surfaces 
whose diameter changes periodically along the cylindri- 
cal axis, in which the diameter function r{z) is given by 
a sinusoidal form. 

For each of these simple surfaces, the Schrodinger equa- 
tion is solved numerically, using a finite difference method 
with appropriate boundary conditions. For instance, a 
Dirichlet or Neumann boundary condition (see Appendix 
IXj) can be used for the case of the catenoid, while the con- 
stant phase increment for Bloch states is used for the case 
of the sinusoidal surface. 



A. Catenoid 

The catenoid is defined in its generalized form by the 
radius function r[z) — a cosh z, in which a (> 0) speci- 
fies the aspect ratio. The surface is minimal if and only 
if a = 1 (see FigQ]), corresponding to the original defini- 
tion of the catenoid. Such surfaces can be thought of as 
simple geometrical models of pore structures common in 
physical systems. Note that the unit of length is taken 
as the characteristic depth of the pore.0 




FIG. 1: The catenoid surface defined by r{z) = a cosh 2 with 
a = 1. The value of a corresponds to the inner radius of the 
pore in unit of the characteristic pore depth along z. 



The variable x defined in Ea. H14|) is given by 

x^ ( \/l + a2 sinh^ z'dz' = -iE{iz, a^), (20) 

^0 



where the function -B((/), k"^) = Jp \/l — siv? (f)d(j) is the 
elliptic integral of the second kind. The potential energy 
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function reads 



a? cosh^ , 



4/ cosh z 4(l + a2sinh zY 



(21) 

where z imphes a function of x. Note that in this expres- 
sion the 'radial' part and the 'longitudinal' part compete 
with each other as a is varied. The radial part, which is 
dominant when a is small, decays outside the pore region 
of unit depth. Without magnetic field, this contribution 
is negative (attractive) for m = while it is positive (re- 
pulsive) for m ^ 0. This is plausible because in a narrow 
part of the catenoid the contribution of kinetic energy 
due to alternation of signs of a wave function around the 
axis becomes costly unless the wave function is effectively 
repelled from that region. This part also inherits some at- 
tractive contribution associated with the curvatures. On 
the other hand, the longitudinal part is dominant when a 
is large, and it decays rapidly when x exceeds 1/a. This 
is purely the effect of the curvature in the longitudinal 
direction and independent of m as well as the magnetic 
field. 

In the case of a = 1, the above expressions reduce to 



f 

Jo 



cosh zdz = sinh z 



and 



1 



1 



1 



(22) 



(23) 



This is the simplest form among the examples we have 
considered apart from a trivial straight tubule. Q Several 
potential energy curves Wmix) for a — I without mag- 
netic field are shown in Fig[21 The minimal value of the 
potential energy for m = 0, Wo(0) = — l/4(l/a^ -|- a^), is 




FIG. 2: Effective potential energy curves Wm{x) for catenoid 
when there is no magnetic field. These four curves corre- 
sponds to m = (bottom), m = ±1, m = ±2, and m = ±3 
(top), respectively. 




FIG. 3: The dependence of the eigenenergies of the bound 
states on the Aharonov-Bohm flux through the pore of the 
catenoid, for (a) a = 0.25, (b) a — 1.0, and (c) a = 2.5. 
Three congruent curves (shifted horizontally) correspond to 
m = —1, 0, and 1. For a = 0.25, there exists the second 
bound state for each m and is shown by a broken line. 



maximal for a = 1, so that the two contributions can be 
said to be balanced for this case. 

We have numerically obtained the energy spectrum 
for the catenoid surfaces with several different values 
of a, where in most cases we have found one or several 
bound states centered around the constricted part of the 
surfaces. 9] The energies of the bound states lie below 
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the continuum > 0. For instance, one bound state is 
found for a = 1.0 without magnetic field and the energy 
eigenvalue has been numerically estimated to e « —0.13. 

For an Aharonov-Bohm type magnetic field, the mag- 
netic field is non-zero only along the rotational axis, so 
that the surface is totally field-free. Hence, any solu- 
tion of Eq. (PJ can be given by a solution of the field free 
equation, multiplied by the phase factor 



exp 



(24) 



which will affect the boundary condition. In the case of 
cylindrical surfaces, this factor is given by 



e<i> 

exp 1 



(25) 



which is independent of z. Note that the function A„i(z) 
in Eq.Q is a constant 



he ' 



(26) 



and the energy spectrum is the superposition of the sub- 
spectra for all the allowed values of A's. This indicates 
that the energy spectrum oscillates as a function of flux 
$ (Fig|3J) and the period of oscillation is given by: 

he 

$0 = — ~ 2.07 X 10""Gauss- TO^. (27) 

e 

In Figl^lwe show the energy eigenvalues of bound states 
as functions of the magnetic flux for a = 0.25, 1.0, and 
2.5. The energy eigenvalue is minimal when the total flux 
is zero and it increases rapidly with the magnetic flux. 
We can also see that the energy spectrum is less sensitive 
for larger a consistent with the fact that the longitudinal 
part of Wmix) is field independent. In the case of uni- 
form magnetic field the oscillation is suppressed for the 
magnetic flux passing through the tubule is no longer 
a constant along the axis direction if the radius is not 
constant. 



B. Sinusoidal tubules 

We next consider the case of infinitely many pores 
along the tubule. This is achieved by taking a radius 
function of the sinusoidal form r(z) = a + bcos{z) with 
a > 6 (> 0), so that the diameter changes periodically 
along the axis. (See Fig01) The relevant variable x is 
given by 

x = [ \J\ + 62 sin^ z'Az' = E{z, -b^). (28) 
Then the effective potential energy curve, 

1 \ 1 1 6^ gQg2 2 



4j {a + bcosz)^ 4 (1 + 62 sin^ z)3 ' 

(29) 



will be periodic in x (see FiglSJ. The solutions of the 
Schrodinger equation for fixed m conform to Bloch's the- 
orem, resulting in the energy bands which depend on 
geometrical parameters of the tubule. The total energy 
spectrum is the superposition of sub-spectra for all m's. 
The consequences of the modulating diameter of the sys- 
tem is evident from potential energy curves. Those for 
the case of constant diameter (straight tubules), these 
functions take merely constant values depending on m, 
posing a constant bias shifts for the energy spectrum. 
In the case of a tubule with modulating diameter, the 
effect of modulation (energy gaps) is prominent in the 
lowest part of the spectrum for each m. The effect of 
modulation can also be seen in the charge density distri- 
bution profiles shown in FigElwith several different Fermi 
energies. Due to the Aharonov-Bohm oscillation of the 
electronic structures, the charge density distribution also 
changes for different values of magnetic flux enclosed in 
the tubule. (FigCJ 




FIG. 4: A sinusoidal surface defined by r{z) 
a — 1.5 and b — 0.5. 



a + 6 cos z with 




FIG. 5: Effective potential energy curves Wm{x) for the si- 
nusoidal surface (a = 1.5, 6 — 1.0) when there is no magnetic 
field. These four curves corresponds to m = (bottom), 
m — ±1, m = ±2, and m = ±3 (top), respectively. 
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FIG. 6: (left) The energy band structures for a sinusoidal 
surface with a = 1.5 and b — 1.0 without magnetic field, in 
which the values of m corresponding to each band is shown in 
parenthesis. Note that for fixed m, energy gaps are the most 
prominent in a lowest part of the spectrum, and it becomes 
dominant for larger m. (right) Charge density distributions 
along the sinusoidal surface with three different Fermi ener- 
gies Ef = 3.0, 3.5, and 4.0, corresponding to the number of 
electrons per period A^e ~ 35.39, 41.49, and 51.61, respec- 
tively. The brighter color corresponds to high concentration 
of the charge density. 



(a) 



Ey =3.0 
1^ =36.9 



(t) 

Er=3.0 
* =0.5*o 
1^=41.02 




¥^4 



FIG. 7: The charge density distribution calculated for Ef — 
3.0 with two different values of magnetic flux, (a) $ = 0.3<l?o 
and (b) $ = 0.5<E>o. 



isfied if a constrained electronic system is realized using 
semiconductor hetero-structures (Af ^ lOOnm) or bent 
carbon graphite sheets[loj| {Xp ^ lOnm). On the other 
hand, the model should be regarded as an excessive over- 
simplification for curved surfaces made of metallic sub- 
stances, in which Xp ~ O.lnm. 



Here we assume the scale of a pore diameter to be of 
the order of a ~ 5nm, which is typically the scale of pore 
structures in mesoporous materials and some of carbon 
nanotubes. Then the magnetic field should be as large 
as B ^ ^ = ^ ^ lO^Gauss = lO^Tesla in order to 
make the total flux passing through the pore to be <i>o. 
This is rather large since the largest static field available 
is about 50 Tesla. The electronic properties could nev- 
ertheless change significantly with more moderate field 
strength if gap openings and/or closings may be induced 
by the magnetic field. It might be interesting to point 
out here that measurements of the Aharonov-Bohm ef- 
fect in carbon nanotubes have been achieved recently. 

[Ulillil 



In our formulation, the initial Schrodinger equation 
reduces to the one-dimensional equations each of which 
corresponds to an angular momentum quantum number 
m along the cylindrical axis. The effective potential en- 
ergy Wm (x) reflects the geometry of the tubule as well as 
the applied field, so that we may actually create a wide 
range of one-dimensional potentials by tuning these fac- 
tors. Since the spectra for different to's are independent 
of each other, a bound state for a m can exist within the 
range of the essential spectrum for another m. We have 
found a simple example with this property in the case of 
a straight tubule with inflated part, r{z) ~ a + b/ (1 + 2:^) 
(a, 5 > 0). It is important to note that the response to 
a change of geometry as well as to an applied field dif- 
fers among different m states, while the Fermi energy 
can cross several partial energy spectra for different rn's. 
One may therefore pursue similar systems for applica- 
tions such as electronic resonators or a selective STM 
probe for different m states. Note that the possibility 
of using carbon nanotubes as STM probes has been re- 
ported, [mil 



V. DISCUSSIONS 

The basic equation (Sec^J assumes hypothetically 
that the surface has zero thickness and the potential en- 
ergy within the surface is flat. One should bear in mind 
that these assumptions pose essential requirement for the 
scale of Fermi wavelength (A^-); that is, it should be suffi- 
ciently larger than the surface thickness and the average 
atomic distance along the surface. This is necessary in 
order for the motion of an electron perpendicular to the 
surface to be strongly quantized and for the flat potential 
approximation to be valid. The condition would be sat- 



If we combine multiple tubules together, more pro- 
found and interesting effects would be expected due to 
the interaction between neighboring tubules. It is also 
interesting to study more complicated surfaces like triply 
periodic curved surfaces Q which are often encountered 
in soft matter systems as well as mesoporous materials. 
In such studies, one may focus on the inherent effects of 
the surface metric, curvatures, symmetry and topology, 
where our present study may also serve as preliminary 
information. We may also pursue the effects of interpar- 
ticle correlation, which should be essential to understand 
real systems. 
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VI. CONCLUSIONS 

We have formulated the problem of quantum parti- 
cle motion constrained on curved tubular surfaces with 
cylindrical symmetry. As the surface has the continuous 
rotational symmetry, the two surface coordinates can be 
separated and the problem is finally reduced to a one- 
dimensional differential equation. Magnetic fields paral- 
lel to the axis is also taken into account in the formula- 
tion. Several examples are studied numerically. 
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APPENDIX A: FINITE DIFFERENCE METHOD 
AND BOUNDARY CONDITIONS 

Using the finite difference method, the original equa- 
tion l(T3|) can be written as 

1 2 1 

— 2yiU - 1)«) + —yi.ja) — 2y((-^' + 1)") 



+W(ja)?/(ja) = ey(ja), (Al) 

where a is the mesh size. Several different boundary 
conditions have been used to obtain numerical solu- 
tions of the one-dimensional Schrodinger equation. First, 
the Dirichlet boundary conditions can be used for the 
case of non-periodic surfaces in which the effective one- 
dimensional potential energy Wmix) behaves asymptot- 
ically as Wm(a;) ^ 0~ (x ^ ±oo). The influence of 
the cutoff ±X of space variable x becomes exponentially 
small for a bound state which is localized in a sufficiently 
smaller region than X, while it will decrease as cx l/AT 
for a scattering state. In both cases, one can estimate 
the energy spectrum fairly well using this boundary con- 
dition with sufficiently large X . If the boundary is far 
enough from the center, so that |W,„(A")| ^ |e|, we can 
improve the precision of the result by taking the Neu- 
mann boundary conditions instead. This means we as- 
sume the form y(X + a) = y(X')e^'^'^ for e — — < 
and y{X -|- a) = yiX^e'^^"- for e = P > 0, where the 
value of e should be determined self-consistently. For the 
case of periodic sinusoidal surfaces, we may use the usual 
assumption of Bloch states, which will lead to a closed 
form of matrix equation. This can be used to obtain the 
band structures. 
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